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ABSTRACT OF THE DISCI nsi jgg 

A system and method for diagnosing the possibility of disease in a 
body part is described. The system includes an electrical data unit for 
measuring electrical data of the body part, the electrical data unit having a 
plurality of N e electrodes. The system also includes a network module for 
representing the body part by a network, the network having external nodes 
and internal nodes connected by current pathways. The system further 
includes an electrical properties module for determining electrical properties of 
the pathways using the measured electrical data, and a diagnosis module for 
utilizing the electrical properties to diagnose the possibility of disease in the 
body part. 
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IRPsis of Disease by. Dgtertninaiion 

Part 



Field of ih f ny —4:~ 

This invention relates to medical diagnosis of disease and specifically 
relates to diagnosis of disease by measuring electrical properties of a part of 
the body. 



Background | mntjnn 

Measurement devices of various kinds have been used to diagnose 
disease. For example, x-ray machines measure tissue x-ray density 
ultrasound machine measure acoustic dens*y, and men™, sensora measure 
differences in (issue heat generation and conduction, which measurements 
have diagnostic value. In addiHon to these devices, .here exist devices that 
can measure etecfrical data, such as vokage or Impedence, of a body part tor 
the purpose of diagnosing disease, su *h as breast cancer. 

Voltage and impedance values,of various types of body tissue ara well 
known threugh /„ Wvo studies on humans or from excised Issue made 
available flowing therapeutic surgtea, procedures . „ is ^ d00umente(( ^ 
a decrease In efeofric- impedance occurs in tissue as « undergoes cancereus 
ohanges. This Ending b consistent over many anlma. species and tissue 
types, including, for example, human breast cancers. 
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However, even though this; relationship between impedance and 
cancer is well known, it is difficult to make straightforward absolute 
comparisons between individuals. For example, if a patient has a smaller 
breast impedance than that of an average female, this does not imply that the 
patient has breast cancer because several fectons could give rise to such a 

smaller value that, are unrelated to cancer. For example, it is known that 

voltage or impedance is related to the percent fat content in an individual. 

Obese women, for instance, would be expected to have breast impedance 

values that are significantly different than average. 

Despite such difficulties, a method that permits comparisons of 
electrical properties for diagnostic purposes has been developed that involves 
homologous body parts, i.e f . body parts that are substantially similar, such as 
a left breast and a right breast. In this method, the impedance of a body part 
of a patient is compared to the impedance of the homologous body part of the 
same patient. 



One technique for screening and diagnosing diseased states within the 
body using electrical impedance is disclosed in U.S. Pat. No. 6,122,544, 
which is incorporated herein by reference. In this patent, data are obtained 
from two anatomically homologous .body regions, One of which may be 
affected by disease. Differences in the electrical properties of the two 
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homologous body parts could signal : disease. One body part, assumed to be 
healthy, is used as a baseline with which to evaluate the homologous body 
part for disease. 

To perform such an evaluation of the body part for disease, a current Is 
injected between an electrode pair and the resultant voltage is measured 
therebetween. This voltage is indicative of the impedance of the body part 
and may have diagnostic value. 

However convenient this method is, obtaining impedance in this 
manner is simplistic. When current is injected, the flow of electricity through 
the body part is expected to take one .or more complex, jagged paths from the 
starting electrode to the final electrode. The conventional approach does not 
account for the internal current pathways that the electricity follows. 

Summary o f the invention 

The present invention describes a system and method for calculating 
the resistive properties of a body part in a way that accounts for the internal 
topology of current pathways. By accounting for this topology, a more reliable 
diagnosis of disease can be achieved. 
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In particular, a system for diagnosing the possibility of disease in a 
body part is described herein. The system includes an electrical data unit for 
measuring electrical data of the body part, the electrical data unit having a 
plurality of N. electrodes. The system further includes a network module for 
representing the body part by a network, the network having external nodes 
corresponding to the location of the plurality of the N e electrodes and internal 
nodes. The internal and external nodes are connected by current pathways. 
The system also includes an electrical properties module for determining 
electrical properties of the pathways using the measured electrical data, and a 
diagnosis module for utilizing the electrical properties to diagnose the 
possibility of disease in the body part. 

For example, the electrical properties module can determine the 
resistance of each of the current pathways, thereby obtaining a conductance 
matrix. The electrical properties module can likewise repeat these steps to 
obtain the conductance matrix of a homologous body part. For instance, the 
conductance matrices of the left and right breast can be obtained and then 
compared by the diagnosis module to diagnose disease. 

Also described herein Is a corresponding method for diagnosing the 
possibility of disease in a body part. The method includes measuring 
electrical data of the body part with a plurality of N e electrodes, and 
representing the body part by a network. The network has external nodes 
corresponding to the location of the plurality of the N e electrodes and internal 
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nodes, the internal and external nodes being connected by current pathways. 
The method also includes determining, electrical properties of the pathways 
using the measured electrical data, and utilizing the electrical properties to 
diagnose the possibility of disease injthe body part. 

Brief description of the drawings , 

Figure 1 is a block diagram of the diagnostic system of the present 
invention; 

Figure 2 is a network of nodes and current pathways, according to one 
embodiment of the present invention; 

Figure 3 is a block diagram of the electrical data unit of the diagnostic 
system of Figure 1; 

Figure 4 is a block diagram of the electrical properties module of the 
diagnostic system of Figure 1 ; 

Figure 5 is a block diagram of the diagnosis module of the diagnostic 
system of Figure 1; 
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Figure 0 is a block digram of the admittance mcdute of the electrica, 
data unit of Figure 3; and 



Figure 7 is a block diagram of the conductance module ofthe electrfoal 
properties module of Figure 4. 

Retailed sjaastiBflgi] ?f«ie inv«n«~> 

Figure 1 shows a diagnostic system 10 for diagnosing the possibility of 
disease in a body part, such as cancer of the breast. System 10 includes an 
electrical data unit 12 having a piMreMy of N. elecbodes ,4, a network 
module 10. an etectrtca. properties mpdute 18, and a dtagnosis module 20. 

EleoMoal data unit 12 measures electrical data ofthe body part using 
the jv, electrodes 14. The V. electees ,4 tncludes electrodes for injecting 
current Into the body part, and electees for measuring resultant vottages In 
one embodiment, the current injection etectrodes and the voltage 
measurement eteCredes coincide. An electrode a n«y containing a plureiHy of 
electrodes similar to one suitable for us. in the present invention is disclosed 
In U.S. Patent No. 0,122,544. 

The electrode array can be connected to e multichannel device that/Is 
ooupted to electrodes 14. A cenbar centre, unit, consisting of a central 
processing unit (CPU, and RAM and ROM memories, seleota, In rapid 



CA 02423445 2003-03-26 



succession, sets of electrodes for current injection and voltage measurement 
One or more analog-to-digital (A/D) converters can convert the analog 
measurements to a digitized form. .Digital data from the A/D converter are 
processed by the CPU, where they undergo real time analyses for error 
checking, routing to a monitor for display of raw or processed data, as well as 
storage in memory for further analysis and output to the monitor and a printer. 

Network module 16 includes hardware and/or software for representing 
y © n©twork has external nodes corresponding to 

the location of the plurality of the N e electrodes 14 and internal nodes. The 
internal and external nodes are connected by current pathways through which 
current flows. The external nodes lie on the perimeter of the network, while 
the internal nodes lie inside. A description of the network appears below in 
connection with Figure 2. 

The electrical properties module 18 includes software and/or hardware 
for determining electrical properties of the current pathways using the 
measured electrical data. For example, the electrical properties module 18 
can use the electrical data obtained by the electrical data unit 12 to find the 
resistances of each of the current pathways in the network. The resistances 
of the current pathways can be arranged as a conductance matrix. 

Diagnosis module 20 utilizes the determined electrical properties of the 
pathways to diagnose the possibility of disease in the body part. For 
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example, diagnosis module 20 canlcompare the conductance matrix of the 
body part to an average matrix obtained from a population group, or to a 
conductance matrix obtained from a homologous body part, as described in 
more detail below. 

Figure 2 shows one possible network .50 of nodes and current 
pathways, according to one embodiment of the present invention. The 
network 50 includes eight external nodes 52, four internal nodes 54 and 
twenty-four pathways 56. The electrical data unit 12 injects current between 
two electrodes and then measures the resultant voltage between two 
electrodes that need not coincide with the two electrodes used to inject the 
current. The electrodes for injecting current and for measuring voltage lie on 
the external nodes 52. 

Electrical properties module 18 then calculates electrical properties, 
such as the conductance of the pathways 56. Network 50 represents the 
body part and the pathways represent possible paths that current travels 
when injected into the body part By using such artwork, a more complete 
description can be obtained of the conductance properties of the body part as 
compared to a network in Which no interna, nodes are included. ,n the 
example of Figure 2, the current pathways are line segments that Intersect 
only at external nodes or internal nodes. 
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Figure 3 shows the electrical; data unit 12 of Figure 1. Electrical data 
unit 12 includes an impedance module 22 and an admittance module 24. The 
electrical data unit 12 can also include an averaging module 25. In circuit 
theory, the impedance. Z is a complex number, whose real part is the 
resistance R and whose imaginary part is the capacitive reactance X c . The 
magnitude of Z is given by V/l, where V and I are the amplitudes of the 
voltage and current, respectively. Both the resistance and reactance can be 
measured, and analyzed in the present invention. 

Impedance module 22 obtains an impedance z v of the body part as 
follows. The electrical data unitinjects a current /, between an ith electrode, 
chosen from the N e electrodes 14, and a reference electrode. For simplicity, 
and without loss of generality, the reference electrode is taken to be the AT e th 
electrode of the N e electrodes 14. The electrical data unit 12 measures a 
resultant voltage V, between a ./th ; electrode chosen from the plurality of 
electrodes 14 and the reference tf,th electrode. The impedance module 22 
calculates Z(f according to z^V i ll ) where.l<:*<;Ar e -l and l<sy 5:^-1. T he 
electrical data unit 12 repeats the step of injecting and the step of measuring 
with varying electrodes to build an N.-lxN e -l impedance matrix Z with 
matrix elements z m given by 

Zm»-VJh (1) 
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where lZmZN.-l and 1& H *N.-1. As described, in more detail below, the 
admittance module 24 calculates art N m -ixN.-i admittance matrix ¥ from 
the impedance matrix Z. 

The averaging module 25 calculates an average admittance matrix, Y . 
In one embodiment, the electrical properties module 18 uses Y to calculate 
the conductance matrix. 

Figure 4 shows the electrical properties module 18 of Figure 1. The 
electrical properties module 18 includes a conductance module 26 for 
determining a conductance matrix, G, for the body part by using the 
admittance matrix Y. The conductance, matrix characterizes the conductance 
properties of the body part Because the presence of some diseases, such as 
cancer, is known to alter the conductance of a body part, the conductance 
matrix possesses considerable diagnostic value. 

Figure 5 shows the diagnosis .module 20 of Figure 1. The diagnosis 
module 20 includes a comparator 28! for comparing the conductance matrix 
for the body part to an average conductance matrix from a population group to 
diagnose the possibility of disease in the body part. A difference between the 
two matrices could Indicate the presence of disease. Alternatively, or in 
addition, the comparator 28 can compare the conductance matrix for the body 
part, such as a left breast, to. a conductance matrix for a homologous body 
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part, the right breast. Again, a difference between the two matrices could 
indicate the presence of disease. 

Figure 6 shows the admittance module 24 of Figure 3. The admittance 
module 24 includes an auxiliary module 30 and an admittance matrix 
calculator 32. 

The auxiliary module 30 calculates an N, -lxN e -1 auxiliary matrix A 
with matrix elements a r The auxiliaiy matrix A is self-adjoint by construction 
and is given by 

» 

A = (Z- l +adjZ" l )/2 ( 2 ) 
where adj denotes the adjoint. 

The admittance matrix calculator 32 calculates the JV e xtf e admittance 
matrix, having matrix elements y v , accoiding to 



a v if !£liSJV,-l andlS/iSJv.-l 



N.-l 



x,Si=N e or y=iv. 



where 

v 

The admittance equation satisfies 



(3) 



*• (4) 



. (6) 
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Where v,is the vector of externa, node potentials whose eiements are the V m 
defined above, and i e is the vector of externa, currents whose elements are 
the ^defined above, in what follows, internal current vectors /, are taken to 
be zero, i.e., current is only injected at the externa, nodes. 

Figure 7 shows the conductance module 26 of Figure 4. The 
conductance module includes a conductance matrix calculator 34. The 
conductance matrix calculator 34. obtains the conductance matrix. G, using 
the admittance matrix by solving the conductance equation 

where, if N . is the number of externa!, nodes and N t is the number of internal 
nodes, then G„ is an ^x^matrix, G el is an N e xN, matrix and G a is an 
AT, xN, matrix defined by 



{gZ.gJ (7) 

The equation G u -g^g^(X) ,fe obtained from the standard nodal 
equation 

by eHminating the vector of internal potentials, v„. from the nodal equation 
using v, =-G?Glv e , and then substituting Yv e = i e . 

Alternatively, the. conductance matrix calculator 34. obtains the 
conductance matrix by solving the conductance equation 



CA 02423445 2003-03-26 



-13- 

where the averaging module 25 calculates an average impedance matrix Y 
by averaging over one or more Impedance matrices, including Y, obtained by 
varying the reference electrode chosen from among the plurality of electrodes 
14. 

The conductance matrix calculator 34 can solve Equation (6) or (8) 
using several methods. Two methods, Newton's method for a system of 
nonlinear equations and a continuation method are now described irfdetail. 

Solution of Co nductance Equation- 

i 

The general form of the conductance equation (Eq. 6 or 8) is 

i 

= -G d CtiG T J-G m =0 (9) 
where f(g) is a function of elements of GG, where g is a vector 
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8= 



8 a 
8* 
8u 
LSuJ 



8* = 



8« = 



8 k = 



8a = 



(9a) 



(A) Newton's Method for a sysfe/n of nonlinear equations 
Let i = l,...,NxN , be a column vector of all possible 

pathways for given Af. Let g 0 be a vector of estimated elements of GG for the 
given network. Vector g has zero for each g {i) that corresponds to certain 

G G{p,q) and nodes p and q are not connected. The Jacobian is jr e -\ = ^ 
Substituting g 0 in (9), /(*,)* o is calculated. 

A new vector g M is calculated from: 



g 



<»+>) _ 



(10) 
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For g 0 close enough to g' (solution of the problem (9)), the iterative method 
described by equations (10) quadratically converges to g\ If go is not a 
good initial guess for g', system (10) may not converge. In that case, the 
continuation method described below is applied. If the number of equations 
is greater than the number of unknowns, a least squares problem is being 
solved using the Gauss-Newton method. In that case, the convergence is not 
quadratic. 

(B) Continuation Method 

For a function /: R" -*r» , a numerical method, sensitive to the initial 
guess, can be used to solve the equation f( g ) = 0 . whereg e R" . if the guess 
is not sufficiently close to the solution g\ the numerical method might not 
converge. To obtain a sufficiently close initial vector $<°> for problem 

=0, an augmented fonction /fetf).*):*"* 1 can be considered 
that is defined as: 



and the following equivalent problem Solved: 

/fe(0),0) = O, 0e[O4] 
where: 



(11a) 
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Xg) :R n ->JT is a function of g, puch that equation /fe)=0 is easy to 

solve, and 0 is a real number. 

The important properties of the new function are: 



/W),0) = 



fig), 0=0 
fig), 6=1 



For 0 = 0, the equation (11a) becomes f(g) = 0, which is easy to solve by 
construction. 

For 0=1, the problem (11a) is identical to the original problem 
/Or) - 0 • For 9 e]0,l[ , the problem.(1 1a) becomes 

a-d)/(g(«))+^(0)) = O 

The following method takes the solution of the problem^ = 0 as the initial 
guess ^)for the close problem a -*)*(*) + 0/Qr(0))=o , where is 0 =0+60 
and S0«l . The underlying assumption is that the solution for the problem 

/Or(0),*)is close enough to g'(0+S0) for the problem hg(e + S8),e+60) so 
that the numerical, algorithm converges to g\e + 39). 
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To summarize, the continuation method includes a) solving the 
problem /(s(0),0) = o for 0=0, b) taking the solution g\0) as the initial 
guess g w (0+8$)for 

Ag(6+se\e+8e) = (L-(e +S6))f(g(6 +se))+(o+se)f(g(e + se)) , 

where 89 is a small step, c) attempting to solve this problem using a Gauss- 
Newton algorithm, d) if the Gauss-Newton algorithm converges, setting 
e = e+Sd, 80 = a-89, where a>l is the acceleration coefficient, and if it 
does not converge, setting <50=<S»/2, ; and e) repeating step b) until 0=1. 

These methods can be used to solve the conductance equation: 

fig) = (G„ -GjGfGl) - G m 
The matrix of finite differences is A/ a/ 2 ... A/J.where 

A(^)=(^-G t e ,c?;V tefc )-g 1Bt , A = i..j\fe. 
From the above equation, 

A/ t -GJ*3?g m -G^Ag^ (14) 



Matrix ag; 1 can be calculated from thp equation =/as follows 
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(15) 



Substituting (15) in (14) leads tp: 



46 = ^ - AG„(G^) + (G^)AG tt (G-^) - (G,^ 1 )^ 



(16) 



The next step is to reorder the matrix products in (16). Each AG„ should 
come as a vector Ag.to the right hand side of the oorresponding matrix 
product. Ag m is a column vector of all columns of AG„ as indicated by (9a). 



By applying the rule 



AXB <=> AAx 



(17) 



Equation (16) becomes 

^ =A S«* + ^ziAga - ^ 3 Ag« (17a) 

where ^ 3 = GefG/7 » , ^ = OJ3 r*^ and M tt and A, are calculated by 
applying (17): 



Equation (17a) can be written as: 



4f t =[Z ~M lk M u -Afj] 



(18) 
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where j k =[/ - Mut M ^ _ M j & the Jacobian matrix corresponding to the 
A* column of F(g< k) ) in equation (10). 



Solving the system 



« 

tat 



*g! n > 



is equivalent to solving equation (19) , 



U 



0 | -A/« | m£> I. -Mj* 
I 5 I : | • 





















Ad* 






J 


.Art*. 







(19) 



To simplify equation (19), the vector of admittances may be written as 

gS*T r - gbmj> (20) 

where 

Sec 

s„ 

8a 
Su 
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Stranch ■» a vector of the current pathways in the selected network whose 
dimension is (Nbx 1), T r is a matrix defined by the network so that equation 
(20) holds: 



The dimension of the matrix T r is (Ne+Ni) 2 xJV6 



otherwise 



Substituting (20) in (19) results in a reduced, system (21): it has the same 
number of equations, but only Nb unknowns and lower rank of reduced 
Jacobian J.: 



(21) 



where 



J = 



I : I I I 

-° ^ i -«s i i o 



-Mi n \ 



and Ag 4rawA is a Vector of current pathways admittance corrections (of size 
Nb). The new admittance vector is calculated from the equation 



(22) 
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Similar to (20), vector can be calculated from g g* as 

A — — * * . ■ ... 

equation (23) 



(23) 

where 7 is a matrix defined by the chosen network so that the 
holds: 



[0 otherwise 
The dimension of the matrix T is ({Ne+NlfxNb). 

There are other methods, besides the method just described, which 
can be used by the conductance matrix calculator 34 for obtaining the 
conductance matrix. For example, the nodal equation introduced above, 
together with Eq. 5 yields 

G m =G„ + G ei v t 

0=^ + 0^ (24) 
where G a sRe(F). , n Equations (24), G m is the measured (i.e. specified) 
conductance matrix and v e does not appear because it is set to the identity 
matrix. In the previous method for obtaining the conductance matrix, 
G le + G uVt was forced to be zero because the internal nodes are not connected 
to the outside, resulting in Eq. 8: 

Equations (24), however, may be solved directly by noting that the W 
could be viewed as a measured >0 m > meaning that it is subject to errors similar 
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to those of G m . This changes the internal nodes from unobservable to 
observable leading to a well-conditioned problem. 

Linearization of Equations (24) gives 

AG a = AG W + AG^v, + G e Av f 

A0 OT = AG- e + AG fl v, + G„Av, < 25 > 
Among the incremental variables, Av, appears and is part of the solution, 
although not of primary interest. The left side of Equations (25) represents an 
extended set of mismatches (A0„ being part). 

Equations (24) may also be solved without linearization by alternatively 
and iteratively solving for the elements in G (with V, being a constant) and V, 
(with the elements of G being held at a fixed value). 

A note about the dimensions of Equations (24) is in order. G is 

nt 

actually obtained by a set of measurements (as many independent 
measurements as there are columns) and for each of these the outside 
voltage is applied to one electrode (plus the base) only. If it is assumed that 
all voltages are given in a per-unit system then the external voltage matrix v, 
may be a unity matrix, hence v e does not appear in the Equations 24. 

Taking note of Equation (8), if v . contains multiple measurements 
arranged so that there is a current only between one electrode and the base 
electrode: 
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where / is the identity matrix. If the system is linear, division by the base 
voltage can be performed and all voltages calculated subsequently are 
expressed with respect to w In this case, the equations 

0=G fe v e + G B v, 
can be re-written as Equations (24). 

Once the conductance matrix calculator 34 finds the conductance 
matrix for the body part, using any of the methods described above, the 
aforementioned steps can be repeated to obtain the conductance matrix of 
the homologous body part. The diagnosis module 20 can then compare the 
conductance matrix for the body part and the conductance matrix for the 
homologous body part by using several comparison methods. For example, 
the norm of the difference of these two matrices can be computed, and if it is 
greater than some threshold, then, further analysis can be performed as this 
difference may signal the presence of disease. 

Several computer systems can be used to implement the method for 
diagnosing disease in a body part. The computer system can include a 
monitor for displaying parts or the whole conductance matrix, or for displaying 
the difference between the conductance matrix for the body part and the 
conductance matrix for the- homologous body part using one of several visual 
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methods. In one embodiment, the method can be implemented on a 2 GHz 
Pentium™ 4 system with 512 MB RAM. 

It should be understood that various modifications could he made to 
the embodiments described and illustrated herein, without departing from the 
present invention, the scope of which is defined in the appended claims. For 
example, although emphasis has been placed on describing a system for 
diagnosing breast cancer, the principles of the present invention, can also be 
advantageously applied to other diseases of other body parts. 
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Claims 

What is claimed is: 

1. A method for diagnosing the possibility of disease in a body part, the 
method comprising 

measuring electrical data of the body part with a plurality of N e 
electrodes; 

representing the body part by a network, said network having external 
nodes corresponding to the location of the plurality of the Jv, electrodes and 
internal nodes, wherein the internal and external nodes are connected by 
current pathways; 

determining electrical properties of the pathways using the measured 
electrical data; and 

utilizing the electrical properties to diagnose the possibility of disease in 
the body part. 

2. The method of claim 1, wherein the step of measuring includes 
measuring an impedance Zy of the body part. 

3. The method of claim 2, wherein the step of measuring an impedance 
Zy of the body part includes 

injecting a current /, between an i* electrode chosen from among the 
Plurality and a reference electrode taken to be the N e electrode of the 
plurality; 
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measuring a resultant voltage V } between a j* electrode chosen from 
among the plurality and the reference; electrode; and 

repeating the steps of injecting and measuring to calculate an 
N e -lxN,-l impedance matrix Z with matrix elements z mn given by 
= K l*n where \<ti<>N e -l andl<S j <ZN e -1. 

4. The method of claim 3, wherein the step of measuring electrical data 
includes calculating an N e ~lxN e -l admittance matrix Y from the impedance 
matrix Z. 

5. The method of claim 4, wherein the step of determining electrical 
properties of the pathways includes determining a conductance matrix, G, for 
the body part by using the admittance matrix. 

6. The method of claim 5, wherein the step of utilizing the electrical 
properties includes 

obtaining an average conductance matrix from a population group; and 
comparing the conductance matrix for the body part and the average 

conductance matrix from a population group to diagnose the possibility of 

disease in the body part. 

7. The method of claim 5. wherein the step of utilizing the electrical 
properties includes 

determining a conductance matrix for a homologous body part; and 
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comparing the conductance matrix for the body part and the 
conductance matrix for the homologous body part to diagnose the possibility 
of disease in the body part. 

8. The method of claim 5, wherein the step of calculating an JV e -lxN e -1 
admittance matrix Y includes 

calculating an N,-lxN,-l auxiliary matrix A with matrix elements a u , 
the auxiliary matrix given by 

calculating the admittance matrix y, with matrix elements y^, 
according to 



x t if i=N e oxj=N e 



n=l 



where 



Xj = 2j a ji- 



9. The method of claim 8, wherein the step of determining a conductance 
matrix includes solving the equation 
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where, if N e is the number of external nodes and N t is the number of internal 
nodes, then G ee is an N e xN e matrix, G ei is an N e xN { matrix and G a is an 
N t xN t matrix defined by 




10. The method of claim 9, wherein solving the equation includes utilizing 
at least one of Newton's method for a system of nonlinear equations and a 
continuation method. 

11. The method of claim 8, wherein the step of measuring electrical data 
further includes calculating an average Impedance matrix F, by averaging 
over one or more impedance. matrices, including. Y 9 obtained by varying the 
reference electrode chosen from among the plurality and analogously 
repeating the above steps. 

12. The method of claim . 11 f wherein the step of determining a 
conductance matrix includes solving the equation 

Gee "GJJ^G^ =Re(F), 
where, if N 9 is the number of external, nodes and N f is the number of internal 
nodes, then G £e is an N.xN, matrix, G ei Is an N e xN t matrix and G a is an 
N i x N t matrix defined by 



> 
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13. The method of claim 12, wherein solving the equation includes utilizing 
at least one of Newton's method for a system of nonlinear equations and a 
continuation method. 

14. The method of claim 1, wherein the current pathways intersect only at 
external nodes or internal nodes. 

15. The method of claim 1, wherein the current pathways are line 
segments. 

1 6. The method of claim 1 , wherein the external nodes lie on a perimeter of 
the network and the internal nodes lieiinside the perimeter. 

17. A system for diagnosing the possibility of disease in a body part, the 
system comprising 

an electrical data unit for measuring electrical data of the body part, 
said electrical data unit having a plurality of N e electrodes; 

a network module for representing the body part by a network, said 
network having external nodes corresponding to the location of the plurality of 
the N t electrodes and internal nodes, wherein the internal and external nodes 
are connected by current pathways; 
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an electrical properties module for determining electrical properties of 
the pathways using the measured electrical data; and 

a diagnosis module for utilizing the electrical properties to diagnose the 
possibility of disease in the body part. 

18. The system of claim 17, wherein the electrical data unit includes an 
impedance module for obtaining an impedance z g of the body part 

1 9. The system of claim 18, wherein, after the electrical data unit a) injects 
a current /, between an I th electrode chosen from among the plurality and a 
reference electrode taken to be the. electrode of the plurality, and b) 
measures a resultant voltage V } between a j ,h electrode chosen from among 
the plurality and the reference electrode, the impedance module 
calculates z 0 according to z 0 =V,/Ij where l£i£N e -l andl^/^-i. 

20. The system of claim 19, wherein, after the electrical data unit repeats 
the step a) of injecting and the step b) of measuring with varying electrodes, 
the impedance module calculates an N,-lxN 9 -i impedance matrix Z with 
matrix elements Zmn given by z ma = V m //„ where l <s i£ N e -1 and 1 <; j < N e -1. 

21 . The system of daim 20, wherein the electrical data unit further includes 
an admittance module for calculating an. jv^-ix^-i admittance matrix Y 
from the impedance matrix Z. 
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22. The system of claim 21 wherein the electrical properties module 
includes a conductance module for determining a conductance matrix, G, for 
the body part by using the admittance matrix. 

23. The system of claim 22, wherein the diagnosis module includes a 
comparator for comparing the conductance matrix for the body part to an 
average conductance matrix from a population group to diagnose the 
possibility of disease in the body part. 

24. The system of claim 22, wherein the diagnosis module includes a 
comparator for comparing the conductance matrix for the body part to a 
conductance matrix for a homologous body part to diagnose the possibility of 
disease in the body part. 

25. The system of claim 22, wherein the admittance module includes 

an auxiliary module for calculating an N e -lxN e -1 auxiliary matrix A 
with matrix elements a n , the auxiliary matrix given by 

A = (z- 1 +adj2T 1 )/2;and 

an admittance matrix calculator for calculating the admittance matrix Y, 
with matrix elements y () , according to 



where 
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a v if l£i£N e -l andl<7<^-l 
x, if i = N e or j —N € 

-f l x n ifi = N t and 7 = ^ 



«=1 



26. The system of claim 25, wherein the conductance module includes a 
conductance matrix calculator for solving the equation 

Gee ~ G ei^it G ei =Re(y) i 

where, If iV c is the number of external nodes and iV, is the number of internal 
nodes, then G ee is an N e xN € matrix, G ei is an N e xN; matrix and G u is an 
N i xN i matrix defined by 

27. The system of claim 26, wherein the conductance matrix calculator 
solves the equation 

by utilizing at least one of Newton's method for a system of nonlinear 
equations and a continuation methods 
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28. The system of claim 25, wherein the electrical data unit further includes 
an averaging module for calculating an average impedance matrix r, by 
averaging over one or more impedance matrices, including Y, obtained by 
varying the reference electrode chosen from among the plurality. 

29. The system of claim 28, wherein the conductance module includes a 
conductance matrix calculator for solving the equation 

where, if N € is the number of external nodes and N t is the number of internal 
nodes, then G^ is an N e xN c matrix, G ei is an N e xN t matrix and G a is an 
N ( xNt matrix defined by 



30. The system of claim 29, wherein the conductance matrix calculator 
solves the equation 



by utilizing at least one of Newtonte method for a system of nonlinear 
equations and a continuation method. 

31. The system of claim 17, wherein the current pathways intersect only at 
external nodes or internal nodes.. 




G eB -G e &Gl=Jte{Y) 
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32. The system of claim 17, wherein the current pathways are line 
segments. 

33. The system of claim 17, wherein the external nodes lie on a perimeter 
of the network and the internal nodes lie inside the perimeter. 
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